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Abstract 

In electroweak baryogenesis, a domain wall between the spontaneously 
broken and unbroken phases acts as a separator of baryon (or lepton) 
number, generating a baryon asymmetry in the universe. If the wall 
is thin relative to plasma mean free paths, one computes baryon cur- 
rent into the broken phase by determining the quantum mechanical 
transmission of plasma components in the potential of the spatially 
changing Higgs VEV. We show that baryon current can also be ob- 
tained using a statistical density operator. This new formulation of 
the problem provides a consistent framework for studying the influ- 
ence of quasiparticle lifetimes on baryon current. We show that when 
the plasma has no self-interactions, familiar results are reproduced. 
When plasma self-interactions are included, the baryon current into 
the broken phase is related to an imaginary time temperature Green's 
function. 
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1 Introduction 



During the electroweak phase transition, baryogenesis can occur in the min- 
imal standard model and its extensions (more Higgs doublets, supersymme- 
try, etc.), provided the phase transition is strongly first order Bubbles of 
broken symmetry (true) vacuum sweep through the unbroken symmetry (un- 
stable) vacuum. The computation of the ratio of baryon density to entropy 
density (measured to be ~ 10^^° 0) proceeds differently if the bubble wall is 
thinQ, ^, ^ or thick[^, |^, ^ relative to the mean free paths of constituents of 
the quark-lepton-boson plasma. In this paper we will be concerned with the 
thin wall case, where one computes the transmission of plasma quasiparti- 
cles past the bubble wall. Because CP and C are violated, the transmission 
of baryon number (or its lepton surrogate) differs from transmission of an- 
tibaryon number. This asymmetry drives the baryon asymmetry of the uni- 
verse, although a number of other physical effects are important in arriving 
at the final result. 

Farrar and Shaposhnikov used this framework to estimate the baryon 
asymmetry produced by the minimal standard model, with CP violation 
originating in the CKM matrix^. They concluded the experimentally ob- 
served number could be generated in this way. This result was criticized by 
Gavela et.ai, who obtained an asymmetry about 10^^'^ relative to that com- 
puted by Farrar and Shaposhnikov |llO[| . These authors state that the main 



reason for the discrepancy is the very short lifetimes of quasiparticles in the 
hot plasma (P ~ 19 GeV for quarks). On the other hand, in the work of 
Gavela et.al, quasiparticle damping is incorporated by giving quasiparticles 
a non-Hermitian effective Lagrangian. This falls short of a complete many- 
body treatment of the self-interacting plasma moving past the Higgs VEV 
(vacuum expectation value) domain wall. 

In this paper we develop an alternate framework for computing currents 
which may make it possible to reduce the number of ad hoc assumptions 
in calculations. We advocate using a statistical operator that specifies the 
thermal system consisting of the hot, self-interacting plasma with the VEV 
domain wall moving through it. In principle, the new framework does not 
require the introduction of quasiparticles at all. It could be used even when 
lifetimes are so short that quasiparticles cease to be useful degrees of free- 
dom. In practice, one must be able to accurately take account of the plasma 
self-interactions, and this usually means working in terms of physical excita- 
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tions, using perturbation theory to account for interactions. Finite quasipar- 
ticle lifetimes and mean free paths are present in the statistical mechanics 
approach as a consequence of self-interactions in the plasma, and when quasi- 
particles decay, baryon current borne by decay products is included. 

Two Lorentz frames are involved in the problem. The "laboratory" frame 
is the rest frame of the plasma, denoted by P. The domain wall moves 
through frame P at speed v. (We take the domain wall to be planar, thereby 
ignoring the curvature of the bubble wall.) The second frame is the rest 
frame of the domain wall, denoted V . In frame V , the plasma streams by 
the stationary VEV wall, from the unbroken into the broken phase. We will 
take the unbroken phase to lie at 2; < 0, and the broken phase at 2; > 0. We 
assume, with references P, ITO], that the latent heat associated with the phase 



transition can be ignored. Then there is no energy source at the bubble wall, 
and we are dealing with a closed system described by a Hamiltonian. Like 
these authors, we also discount turbulence so that the system is stationary 
in frame V . Taken together, these observations imply that the plasma is 
described by observers in \^ by a time-independent statistical operator p. 
Because the thermal bath in contact with the gas is at rest in frame P, the 
density matrix is 

p = exp[-/3irp], (1.1) 

where Kp = Hp — fiN, and Hp is the Hamiltonian in P. What forces us to 
deal with two frames is that while Kp determines p, Kp is time-independent 
only in frame V. [Recall that to recast the computation of the baryon current 
as a conventional problem in equilibrium statistical mechanics, p must be 
time-independent.] 

In frame V there is a conserved baryon current J^(x, t) whose form and 
time-dependence are determined by Hy, the Hamiltonian of the system in 
V.f\ The baryon current past the VEV wall, measured in V, is 

This expression appears simpler than it is because one Hamiltonian appears 
in p, a second Hamiltonian controls time evolution, and the two do not 
commute. As a result, J generally depends on z and t. However, as t ^ cxd. 



is used when we wish to distinguish an operator from a related c-number O. 
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transients decay and the current approaches a hmit that is constant in time 
and space. In the case of a plasma without self-interactions, the late time 
current is what is expected on the basis of simple physical arguments. The 
resulting formula (Eq. |3.16| ) has been used in previous work, either implicitly 
or explicitly. 

The matter content of the plasma in the epoch of baryogenesis depends 
on whether the standard model or one of its extensions is under study. Even 
were we to choose a specific case, the inclusion of many species would tend 
to obscure the gist of our work. In a similar vein, implementation of real- 
istic baryon-nonconserving processes would introduce complexity and detail 
that are undesirable. What is illuminating is to study a simple system for 
which the statistical matrix approach is developed in a way that can be 
generalized to other plasmas. In this paper we assume a single spin zero 
baryon-antibaryon species is present in the plasma. The Higgs VEV intro- 
duces a z-dependent mass m{z) for these bosons. We assume that the VEV 
jumps to its broken symmetry value only in the region < z < L. That is, 
in this interval, mlz) ~ M, while elsewhere m{z) ~ m, and m < M. Having 
m{+oo) = m(— cxd) has the advantage of simplifying the asymptotic form 
of basis functions that appear extensively in the calculations, and it is for 
this reason that the width of the broken symmetry phase is finite. Since L 
can be as large as we wish, its finitude has no physical significance because 
tunneling through the broken symmetry phase can be made as small as we 
wish by increasing L. 

If nothing further were done, we would find J = because while there 
would be a baryon current past the domain wall, there would be an equal 
and canceling antibaryon current. We produce a net baryon current by in- 
troducing a chemical potential that enhances the baryon density relative to 
the antibaryon density. It plays the role of baryon non- conserving processes 
in our model. In order for the Bose distribution functions to be finite, the 
chemical potential must lie in the range —m < ^ < m. This restriction is 
what requires us to choose m 7^ in the "unbroken symmetry" phase in our 
model. 

In Section 2 we consider the diagonalization of Kp, Hp, Ky and Hy-, in 
the case where the plasma has no self-interactions. This allows us to compute 
the baryon current of such a plasma in Section 3. In the limit t — 00, the 
baryon current we find is the expected one. In Section 4 we consider the 
case of a plasma with self-interactions. The objective there is to express 
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the current in terms of an imaginary time temperature Green's function in 
equilibrium many-body theory. It is these Green's functions that have simple 
expansions in perturbation theory, and to which the formalism and methods 
of many-body theory apply. In Section 4 a number of complications come 
together: the difference between Hp and Hy, the need to take t oo, and 
the requirement of relating real time and imaginary time temperature Green's 
functions. We succeed in reconciling these strands only by working to first 
order in the plasma speed v. (In Section 3 our results for a plasma without 
self-interactions hold for < f < 1.) We argue that the baryon current J{v) 
should be odd in v, so this restriction is less stringent than it might seem. 
Extension of our results to 0{v^) is conceivable, but would be complicated. 
At the end of Section 4 we have a specific imaginary time temperature Green's 
function which gives the baryon current to 0{v). The effect of quasi-particle 
decay on baryon current is now accounted for by assessing the effect of plasma 
self-interactions on this Green's function. In particular, contributions to the 
boson propagator where a boson decays on one side of the VEV wall and 
reassembles on the other correspond to contributions to the baryon current 
from quasi-particle decay products. 



2 Free Fields 

To show that a statistical formulation of the calculation of baryon current 
works, we first consider a plasma that has no self-interactions, but full inter- 
action with the VEV wall. As a preparation, we formulate and diagonalize 
the four operators KpQ, Hpo, Kyo, and Hyo- Hyo, the free Hamiltonian in 
frame V, differs from the conventional free boson Hamiltonian only through 
its spatially varying mass term: 



vo 



drx 



(2.1) 



HvQ generates ^-displacements; we also need the momentum operator Pv,z-i 
which generates ^-displacements: 



V,3 



(fx 



dz 



dz 



(2.2) 
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Then for any Heisenberg operator O, 

do do 

— = -i[0,Hvo]; -^ = ~i[0,PvA- (2-3) 
Observers in frame P attribute primed coordinates to events, where0 

t' = -f{t-vz)] z' = -f{z-vt). (2.4) 

Thus, 

dO dO dO . 

_ = ^_ + = [o, ^Hvo + IvPvA . (2.5) 

Since the generator of t'-displacements is HpQ, we have 

Hpo = -fHvQ + jvPv,3- (2.6) 

(A similar formula holds for a plasma with self-interactions; it will be used in 
Section 4.) Finally, the grand canonical Hamiltonians are Kyo = Hyo — l^N 
and KpQ = HpQ — nN, where the baryon minus antibaryon number operator 
is 



N = i (Px 



(2.7) 



KpQ is the key to diagonalizing these operators, since once Kp^ is in 
diagonal form, the remaining operators follow by setting f = 0, or yU = 0, or 
both. The construction begins by considering the Heisenberg equations of 
motion that are generated by KpQ. 

-^ = -#-#-^0- (2.9) 

7 at oz 7 

This equation implies that the current is conserved, where 




= _iy0tv^0, (2.10) 



is the relative speed between P and V. We use /3 for 1/kT. 
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and satisfies Eq. |2.^ . (When we replace the free Hamiltonian Hpo by the full 
Hamiltonian Hp, Eq. p.8| acquires additional interaction terms and countert- 
erms on the right. is still conserved in the theory with self- interactions.) 
Setting f = yU = 0, the conserved baryon current in frame V and Eq. is 



Jo 



- d 



(2.11) 



We obtain basis states for the diagonalization of KpQ by finding classical 
solutions of Eq. plSj If we write (p = 0exp[i/xt], the chemical potential drops 
out, and d> satisfies 



2'~fv 



7 m (z) 



(2.12) 



dt^ ' dz^ ' dtdz 

By separation of variables, we have solutions = exp(— 'ii?t)P(k, x) where 



exp[ikx-xx + z7^^2;] 
^(k,x) = Zp[k^,z), 



E = Vk^ + m'^, and 

d'^Zp{k3,z) 



dz^ 



[2E(27r)3]i/2 



+ YiH + m'- m\z)]Zp{ks, z) = 0. 



(2.13) 



(2.14) 



The current defined in Eq. p.lO| continues to be conserved when is re- 
placed by one classical solution, and 0''^ by the complex conjugate of a second 
solution. Consider the integral of the resulting charge density. Introducing 
useful notation, we write the total charge as 



(0(k),0(kO)=e^(^-^')*(P(k),P(k')), 



(2.15) 



(P(k),P(k')) = j rf^xP*(k,x) 



E + E' — i'fv-^ 
oz 



Att^EE' 



dze'^''^^'-^^UZ*p{k3,z) 



-f{E + E') - iv 



P(k',x 

A 

dz 



Zp{k'^,z) 



However, the total charge should be time independent, implying that 
(P(k), P(k')) vanishes for E ^ E' . This makes (P(k), P(k')) an appropriate 
inner product for the time-independent solutions. 
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Orthogonality can be verified directly. It follows from Eq. p.l4| that 



d_ 

dz 



(2.16) 



3, z) 



-f{E + E') - iv 



d_ 

dz 



Integrating this, the left side vanishes, so the inner product is zero unless 
E = E' or k3 = ±k'^. 

When ks = ±^3, there is a delta function contribution to the inner prod- 
uct generated by the integral over z. Its strength is determined by the asymp- 
totic behavior of Zp, which we now specify. At large \z\, solutions of Eq. |2.14 
are plane waves, and we adopt scattering boundary conditions. For k^ > 0, 

2 (k z) = I ^^P(^^^32) + rpiks) exp{-i-fksz), {z « 0), ,^ 

^' ' \tp{k3)exp{i-fksz), {z»L).' ^' ^ 



and for ks < 0, 
Zp{k3,z) 



r tp{k3) exp{i-fk3z), {z « 0), ,^ 

\exp{i'yk3z) +rp{k3)exp{-i'yk3z), {z » L). 



These scattering amplitudes satisfy "unitarity" relations that we use repeat- 
edly. Eq. |2.14| implies that the expressions 



Zp{k3,z) 



d_ 

dz 



Zp{k'^,z), Zp{k3,z) 



d_ 

dz 



Zp(-k3, z) 



are independent of z when k'3 = ±k3. Evaluating the expressions at large 
positive and negative z, we obtain the unitarity relations 



\rp{k3)\' + \tp{k3)\' = l, 
rUh)tpi-k3)+t*p{k3)rp{-k3) = 0, 
tp{-k3) = tp{k3). 



(2.19) 



Using Eq. p.lS , the complete orthogonality relation is 

(P(k),P(k')) = 75(k-k'). 



(2.20) 
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Solutions studied so far have E > 0. Related negative frequency solutions 
are the complex conjugates = exp(+z£'t)P*(k, x). The inner product of a 
negative energy solution with a positive energy solution is 



((P*(k),P(k')) = 1 d^xP{k) 



E' - E- i-fv 



d_ 

dz 



(2.21) 



75(k± + k'j 



7(E' - E)-iv 



d_ 

dz 



Zp%,z) 



It can be shown that the integral over z always vanishes. The complete set 
of orthonormality relations is 

(P(k),P(k')) = -(P*(k'),P*(k)) =75(k-k'), (P*(k),P(k')) = 0. 

(2.22) 

The quantum fields can be expanded in terms of these solutions. 



(2.23) 



0(x,t) = e'^*y d'^k [e'^^*a(k)P(k,x) + e'^V(k)P*(k,x) 
Using Eq. we find 

/■ _ r _ / f)\ 

Pfk,x) 



(2.24) 



7rt(x,t) = -le""' j d^k e-'^*a(k) - 

-e^^V(k) (^^ + iv-^^ P*(k,x) 

Now set t = to obtain Schrodinger picture fields. Using orthonormality, we 
can project out the operator coefficients 



a(k) 



d^^X 



(x) + iv-^) P*(k, x) + Z7rt(x)P*(k, x) 



(2.25) 



b\k) = j d^x 0(x) |^^-it;^^P(k,x)-i7rt(x)P(k,x) . 

Canonical commutation relations for the fields imply that these operators all 
commute, except 



(k),at(k')l = [6(k),6t(k')l =(5(k-k'). 



(2.26) 
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These operators destroy baryons (a) and antibaryons (b). 

The reason our construction leads to a basis that diagonahzes Kpo is that, 
by Eqs. f^.23| and ^.24| , the time dependence of the destruction and creation 
operators is sinusoidal. This implies 



[a(k), Kpo] = {E- /i)a(k), [6(k), Kp,] = {E + ^)6(k). 



(2.27) 



For these equations to hold we must have, aside from a normal ordering 
constant that is irrelevant in this problem. 



K 



PO 



\e - fi) a^(k)a(k) + {E + fi) b\k)b{k) 



(2.28) 



This equation may be derived directly by substituting the expansions for the 



Schrodinger picture operators into Eqs. and |277 . 

When V = 0, KpQ becomes Kvq, and the basis functions P(k, x) are 
replaced by 

exp[zk_L ■ x_|_] 



y(k,x) 



\2E{2 



TT 



(2.29) 



+ [A;^ + - m^{z)]Zv{k3, z) = 0. 



The fields may also be expanded in this basis: 



(x,t) 



7rt(x,t) 



d^k 



e-'^'A{k)V{k, x) + e'^'B\k)V*{k, x) 



(2.30) 



d^kE 



e-'^'A{k)V(k,^) - e'^'B\k)V*{k,^) 



The new operators A and B also have the commutation relations to destroy 
baryons and antibaryons. The operators Kyo and Hyo are diagonal when 
expressed in terms of these operators: 



Kvo= I d''k[{E-fi)A^(k)A(k) + {E + fi)B^{k)B(k) 
Hvo= I d^kE[A\k)A{k) + B\k)B{k) 



(2.31) 
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3 Current of Free Baryons 



For free baryons, the obstacle to the evaluation of Eq. |1.2| is that p is expressed 
in terms of frame P operators a, . . ., while the Heisenberg picture operator 
J3(x, t) is expressed in terms of frame V operators A, . . .: 



J3(x,t) 



,i(£-E')t^t(k)A(k')l^*(k,x) V(k',x) 



+e^(^+^')*At(k)5t(k')F*(k,x) ^-gj V*{k',^ 
+e^(-^-^')*5(k)A(k')V^(k,x) (£) V{k- ■ 
+e^(-^+^')*5(k)5t(k')K(k, X 



(3.1) 



l^*(k',x) 



Because E > 0, Wk^, oscillations must suppress the contributions of the two 
middle terms as t — oo. To evaluate the contributions of the remaining 
terms, we must write the A, ... in terms of the a, . . .. Begin with the projec- 
tions 



^(k) 
5t(k) 



(fx 
(fx 



[x)EV^*(k, x) + i7r^(x)V*(k, X 
[x)El^(k,x) -Z7r^(x)F(k,x 



(3.2) 



Substitute expansions |2.23| and |2.24| , evaluated at t = 0: 

/oo 
dk'^ [a(k^ + zk'^)Ri{k^,k3,k'3) 
-oo 

+6^(-kx + zk'^)R2ik±, A;3, /c^)] , 

/oo 
dk'^ [a(-kx + zk'^)Rl{k^, k3, k'^) 
-oo 

+b^{k^ + zk'^)Rl{k^,k3,k'^ 



(3.3) 



where 

Ri{k^,k3,k3') 

R2ik^,k3,k3') 



.oo dze'^"^'^ / 

/■oo ^^p-i'yvE'z ( - ^ 



(3.4) 

Zp{k'^, z). 
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One of the traces we need for the baryon current is 



Tr 



pA^{}^)A{V)\= Too dq,dq'. 



(3.5) 



X 



+R*,{k^, ks, q3)R2ik'^, k'„ q!,)Tr [p6(-k^ + %3)&^(-kl + zq',)] } 



The traces in Eq. |3.5| are standard [11], leading to the expression 



Tr 



pA^{k)A{k') 



Tr [p] 



SA{k±, ks, k'^ 



dq3 



= 5(kx-ki)5^(A;x,fc3,4), 
Rlik±, ks, q3)Ri{k±, k'^, q3) 



(3.6) 



+ R*2ik±, k3, q3)R2{k±, /cg, gs) 



1 



In this equation, E = yk^ + q'^ + m?. There is a similar formula for 
Tr[pB(k)B'^ {k!)]. Sb differs from Sa by the interchange Ri 
The baryon current at late time is 



R, 



J{z,t) 



(Pk^dk^dk'^ 



e'(^-^')*5^(fc^, fcs, fc^) - e-(^-^')*^B(fc^, ^3, k'. 



(3.7) 



To extract the limiting current, consider the contribution of the term pro- 
portional to 5*^. First integrate over k^, for fixed values of k^ and fcg. The 
k^ integration is initially along the real axis of the complex k^ plane. Deform 
the contour to C, which runs from k^ = —oo slightly below the real axis 
(in the third quadrant), crosses the real axis at k^ = 0, and continues to 
k^ = +00 slightly above the real axis (in the first quadrant). Except at the 
origin, contributions on C are exponentially damped as t — > oo, and at the 
origin contributions are suppressed by oscillation. The contribution to the 
limiting value of J{z, t) is therefore zero unless singularities are encountered 
when the contour is changed from the real axis to C. A time-independent 
contribution arises if 5*^ has a pole bordering the real axis, a.t k^ = \k'^\ + irj, 



or at k^ = —\k'^\ — irj, since at k^ 



ibfcg, we have E 



E'. In fact Sa as 
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defined by Eq. p.6| has such poles. Note that when ks = it/cs, Eq. |2.14| imphes 
that J is also independent of z. 

By Eq. |3.(j| , poles in Sa_ arise where poles of Ri pinch the ^3 integration 
contour. The Ri have poles adjacent to the real axis owing to divergence of 
the integrals at z = ±00 in Eq. |3^ . The positions and residues of these poles 
are determined by the asymptotic forms of Zy and Zp, and therefore depend 
on the profile of the VEV wall only through the scattering amplitudes. 

The asymptotic behaviors change with the sign of the real parts of k^, k'^, 
and qs, and many contributions must be added to determine the residues of 
the poles of Sa at k^ = ±^3. Here we consider only two of the contributions 
to illustrate what is involved. 

When ks > and ^3 < 0, the poles of Ri arising at 2; = —00 are the same 
as those of 

Ri{k^,ks,qs) ~ / ^^y>'s^ + ry{ks)e-''^^' 



X (e, + ^E, - IV j-^ tp{q,)e'^'^^^ (3.8) 
tp{q3){Ek + 'jEq + 'jvq3) 



+ 



Am^EkEq 



-fqs + -fvEg -ks-ir] 



7^3 + jvEg + ks-ir] 



The first denominator vanishes at ^3 = 7/^3 — 'jvEk, and the second at ^3 
—7^3 — 'jvEk- We rewrite the terms to exhibit poles in the ^3 plane. 



Riik±,k3,q3) 



tpilks - -IV Ek) hEk - -fvks 



2iTc{q3--fk3 + -fvEk-ir])\l Eu 



tp{.-lh - lvEk)r\,{k3) hEk + ivk. 



3 



2m {q^ + 7/C37 + vEk - ir]) \ Ek 



(3.9) 



Substituting either of these contributions into Eq. ^]6| produces a pole at 
k^ = k'^. The two contributions are 
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SA{k±,k3,k'^)^ (3.10) 



2niE'^{^k3 - -fvE^ - -fk'^ + -fvE'^ - ir])[exp P{-fE'^ - -fvk'^ - //) - 1] 



^ \tp{lk', + ^vE'Mrv{k',)\\jE', + 7.;^^) 

2mE'^{-^k3 - ^vEk + ik'^ + ^vE' - iri)[exp P(-fE'^ + ^vk'^ - A*) - 1] ' 

The step function factor in the first term imposes ^3 < so that the assumed 
asymptotic behavior holds at the pinch. Rewrite these terms to exhibit the 
poles in the k^ plane. 

SA{k^,k3,k'^) - (3.11) 
2m{k3 -k'^- ir/)[exp/3(7E^ - 7^^ - jj) - I] 



\tp{lk', + ^vE',)\yy{k',)\'' 



2m{k2. -k'2 + i7])[expp{'yE'^ + 'yvk'^ - /x) - 1] ' 

The first of these poles is in the first quadrant, and contributes when the 
contour is deformed to C, but the second pole is in the fourth quadrant and 
is irrelevant. 

When all contributions are added, and the unitarity relations invoked, 
the relevant singularity of Sa is 

/> > sgn(fc^) 

^A[K^, Ks, K,, ~ 2^.^^^ _ ^, _ ^^3g^(^,)][g^p^(^^, _ - ;,) - 1] • 

(3.12) 

There is no relevant pole at k^ — —k'^. The term Sa makes the time- 
independent contribution 
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Ja 



[exp P{^E - ^vks - //) - 1] 



(Pkksltviks 



1 



(27r)3E [exp P{-fE - ^vk^ - /x) - 1] ' ^^'"^^^ 

In the term proportional to Sb, deform the fcs- integration from the real 
axis to contour C", which lies in the second and fourth quadrants of the 
complex /ca-plane. There is then a time-independent contribution from the 
following pole in jS'b: 



SB{k±,ks,k'^) 



(3.14) 



sgn(A;0 



27ri[k3 -k's + ivsgn{k's)] { [exp P{-fE'f^ - -fvk'^ + fi) - I] 
The addition to the baryon current is 

1 



+ 1 



Jf 



-I 



(Pkks\tv{ks)\^ j 

(27r)3£; \ [exp Pi^E - ^vk^ + //)-!] 



+ 1 



(3.15) 



The extra term in the brace is present because we did not normal order 
the current operator. No matter: The resulting integrand is odd in k^ and 
vanishes. Altogether, the baryon current in the limit t — > oo is 



J 



-J 



d^kks\tv{ks)\' 



{2ttYE \ exp I3{^E - -fvk^ -pi)-! 



exp f3{'jE - 'jvk^ + At) - 1 j ■ ^^'^^^ 

This result is expected. The Bose distributions are those of a moving 
free gas; they give the densities of baryons and antibaryons, which differ for 
/X 7^ 0. The contributions to the net current density are opposite for baryons 
and antibaryons. They are given by particle density multiplied by velocity, 
ks/E, and the quantum mechanical probability of transmission past the VEV 
wall. 
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4 Plasmas with Self-Interactions 



The result of the last Section shows that the approach based on statistical 
mechanics works for a plasma without self-interactions, for any velocity of 
the plasma relative to the VEV wall. The methods used there relied on the 
absence of self-interactions. On the other hand, the physical question of the 
influence of quasiparticle lifetime on the baryon current can be addressed only 
when we take self-interactions into account in a systematic way. This can 
be done by reformulating the problem as a standard problem in many-body 
theory. In particular, we want to reduce the problem to the computation of 
a temperature Green's function. 

In this Section, we so recast the problem, but only to first order in v. 
Note that the current of Eq. is odd in v: J{—v) = —J{v). We expect 



this reflection property to hold in a self-interacting theory also, when the 
VEV wall profile in the vicinity of 2 = L is the mirror image of the profile in 
the vicinity of 2; = 0. Then the correction to our 0{v) result is O(f^), and 
the 0(y) result is applies over a useful range of v. 

Write the statistical matrix of the self-interacting plasma in the form 

p = e-^^^ = e-'^^^^U{(3y, U{t) = g-T^v-g-rKp^ ^^^^^ 
Introduce the boost operator, and its "imaginary time Heisenberg" picture: 
B = -fKv - Kp; B{t) = e"^^^ Be'^^^^ . (4.2) 

Then 



dU{T) 

dr 



B{t)U{t)- (4.3) 



00 1 T- 

U{r) = E^/ dn---dTnTr[B{n)---B{rn)], 

where in the product, larger r operators stand to the left of smaller r oper- 
ators. 'jKv and Kp have the same self-interactions, so 



B = -'yv J d^x 



— — TT' + Ti — 

oz oz 



(4.4) 
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The power series for U is therefore equivalent to an expansion in powers of 
V. Through first order in v (but all orders in self-interactions) 

JTr[p\ = Tr[e-^^^ J3(x, t)] + dTTr[e~^^'' BMM^, t)] + . . . . (4.5) 

Jo 

The first term on the right is the baryon current of a plasma at rest in frame 
V; it is zero. The second term is the first order baryon current J(l,x, t): 

J(l,x,t)Tr[p] = f dTTr\e-^^^B{T).h{-KM. (4.6) 
Jo 

The Lehmann representation for J(l,x, t) is |12 
J(l,x,t)Tr[p] = Y.j dEmdEnc'^''--''-^' 



(4.7) 



X / danidan{rn\Xi{yi)\n){n\B\m). 



The sum is over eigenstates of Hy and N . It is partially discrete and partially 
an integral over dEda; the state energy dE is exhibited. We expect the inner, 
partially integrated, expression to include a contribution proportional to an 
energy delta function: 

J damdan{m\J3{x)\n){n\B\m) = 5{E^ - En){m\J\n) + ... (4.8) 

The delta function contribution arises because operator B is an integrated 
density. It is present for a free plasma, as we will show below, and it is there 
in each order in the perturbative expansion in the plasma self-interaction. 
Furthermore, current conservation requires {m\J\n) to be independent oi z, 
because Em = En- Between eigenstates of ify, current conservation states 

dJ 

= (^1^ + V ■ J|n) = {m\i{Em - En) Jo + V • J|n). (4.9) 
When Em = En, this becomes 

d 

= V-(m|J|n) = — (m|J3|n). (4.10) 
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Only the energy- conserving contribution survives in the hmit t oo, leaving 
the finite ^-independent current 



J{l)Tr[p\ = J(l,x,t = oo)Tr[p] = / dE^{m\J\n) 

m,n 



(4.11) 



En — Ern 



We can relate J(l) to the imaginary time temperature Green's function 
J3(x,r) = e"^^J3(x)e-"^^. (4.12) 



This Green's function has familiar properties: It is a function of ri — r2, and 
is periodic in this variable with period f3. The Fourier transform ^(x, cUp), 
with ujp = 2ttp/P, is 



(4.13) 



This Fourier transform is the amplitude that has a simple expansion in Feyn- 
man diagrams in many-body perturbation theory. Its Lehmann representa- 
tion is 



g-/3(£;„-Ar™At) _ ^-l3{Em-NmH) 



E — E 



X / damdcrn{m\j3{x)\n){n\B\m). 



(4.14) 



We see immediately that — ^(x, tUp = 0) is identical with J(l,x, t = 0). 
However, Eq. [4.14| has contributions from states with Em 7^ En that make it 
differ from J(l), Eq. [4.11| . Fortunately, there is a simple way to eliminate the 
contributions from unwanted states and recover J(l). We see from Eq. ^]6| 
that states with E^ = En make contributions to Q{'x,ujp = 0) that are 
independent of z, while the unwanted states make contributions that depend 
on z. There is no difficulty in extracting the 2;-independent contribution; it 
pops out when the integration in Eq. is carried out. Thus the relation 
between J(l) and the temperature Green's function is 



J(l) 



-^(x,a;p = 0)|^_i„d. 



(4.15) 
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Note that no analytic continuation in uj is required in this problem, in 
contrast to the situation when one studies real-time transport phenomena at 
finite temperature. 

We are primarily interested in ^(x, tUp = 0) when the plasma has self- 
interactions. However, it is straightforward to compute it for a free plasma. 
It is useful to do so, because the emergence of the energy conserving delta 
function is illustrated, and the correctness of the result can be confirmed. 
We find 



-^(x,u;p = 0) 



d'^k±dk^dk'^ 



{E - E')[e'^(E-^') - l][e^(^'-A') - 1] 
~ {E - E')[ef^(^+>') - l][e/3(^'+/^) - 1] 



(4.16) 



roo 

X / dz' 

-oo 



.„dZy{k'^,z') , I . , dZv{k2„z') 
dz' ^^(^3, z)+tE Zy{k^, z ) — 

A{2'kYEE'{E + E')[e/5(^-M) - i][eP(E'+iA _ i] 



X 



dz 





POO 




/ dz' 

J —oo 



+ 



Zv{ks,z) I I Zvik'^,z) 



dz' 



dZv{k'^,z') 

dz' z ) 

, , , , dZv{k3,z') 
-iEZv{k^,z) — 

iE^-^y^^Zl.{k,,z') 



+^E'Zm,z')^^^y^ 



First note that this Green's function does not depend on x^. The reason is 
that the transverse integrations in Eq. [4.4| set k_L = k'j_. However, k^ ^ k'^ be- 
cause of the VEV wall, and the Green's function certainly depends on z. The 
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integrations over z' come from Eq. [4.4| . They do not impose momentum or 
energy conservations generally, but the asymptotic tails of these integrals do 
produce energy conserving delta functions in addition to other contributions. 
Such delta functions lead to discrete position-independent contributions to 
^(x, cijp = 0). The delta function terms occur in just one of the integrals: 



dz' 



(4.17) 



k',) + r*yi-k,)tyiks)5ik, + k',] 



+ .... 



Retaining just the delta functions, the z-independent contribution to 
^(x, tUp = 0) implies 



J(l) = -^(X,u;p = 0)|^_ind 
d'kkl\ty{k,)\' I 



(4.18) 



1]2 [eP{E+t,) _ 1]2 



This agrees with the 0{v) baryon current obtained by expanding Eq. p. 16 



Physical questions about the baryon current are thus converted into ques- 
tions about a temperature Green's function. We have experience and tech- 
niques for the study of such objects. 
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